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Abstract — The concept of the analytic signal is extended from 
the case of a real signal with a complex analytic signal to 
a complex signal with a hypercomplex analytic signal (which 
we call a hyperanalytic signal) The hyperanalytic signal may 
be interpreted as an ordered pair of complex signals or as a 
quaternion signal. The hyperanalytic signal contains a complex 
orthogonal signal and we show how to obtain this by three 
methods: a pair of classical Hilbert transforms; a complex 
Fourier transform; and a quaternion Fourier transform. It is 
shown how to derive from the hyperanalytic signal a complex 
envelope and phase using a polar quaternion representation 
previously introduced by the authors. The complex modulation 
of a real sinusoidal carrier is shown to generalize the modulation 
properties of the classical analytic signal. The paper extends 
the ideas of properness to deterministic complex signals using 
the hyperanalytic signal. A signal example is presented, with its 
orthogonal signal, and its complex envelope and phase. 



I. Introduction 

The analytic signal has been known since 1948 from the 
work of Ville |[T| and Gabor |2|. It can be simply described, 
even though its theoretical ramifications are deep. Its use in 
non-stationary signal analysis is routine and it has been used in 
numerous applications. Simply put, given a real-valued signal 
f{t), its analytic signal a{t) is a complex signal with real 
part equal to f{t) and imaginary part orthogonal to f{t). 
The imaginary part is sometimes known as the quadrature 
signal - in the case where f{t) is a sinusoid, the imaginary 
part of the analytic signal is in quadrature, that is with a 
phase difference of —7t/2. The orthogonal signal is related 
to /(t) by the Hilbert transform |4,|. The analytic signal 
has the interesting property that its modulus \a{t)\ is an 
envelope of the signal f{t). The envelope is also known as 
the instantaneous amplitude. Thus if f{t) is an amplitude- 
modulated sinusoid, the envelope |a(t)|, subject to certain 
conditions, is the modulating signal. The argument of the 
analytic signal, Za{t) is known as the instantaneous phase. 
The analytic signal has a third interesting property: it has 
a one-sided Fourier transform. Thus a simple algorithm for 
constructing the analytic signal (algebraically or numerically) 
is to compute the Fourier transform of f{t), multiply the 
Fourier transform by a unit step which is zero for negative 
frequencies, and then construct the inverse Fourier transform. 

In this paper we extend the concept of the analytic signal 
from the case of a real signal f{t) with a complex analytic 
signal a{t), to a complex signal z{t) with a hypercomplex 
analytic signal h{t), which we call the hyperanalytic signal. 
Just as the classical complex analytic signal contains both the 
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original real signal (in the real part) and a real orthogonal 
signal (in the imaginary part), a hyperanalytic signal contains 
two complex signals: the original signal and an orthogonal 
signal. We have previously published partial results on this 
topic |5|-|7J. In the present paper we develop a clear idea of 
how to generalise the classic case of amplitude modulation to 
a complex signal, and show for the first time that this leads 
to a correctly extended analytic signal concept in which the 
(complex) envelope and phase have clear interpretations. We 
show how the orthogonal signal can be constructed using three 
methods, yielding the same result: 

1) a pair of classical Hilbert transforms operating indepen- 
dently on the real and imaginary parts of z{t), 

2) a complex Fourier transform pair operating on z{t) as a 
whole, with positive and negative frequency coefficients 
rotated by ±7r/2 using i times a signum function, 

3) a quaternion Fourier transform with suppression of neg- 
ative frequencies in the Fourier domain. 

All three of these methods are based in one way or another 
on single-sided spectra, but it is only in the third (quaternion) 
case that a single one-sided spectrum is involved. 

The construction of an orthogonal signal alone would not 
constitute a full generalisation of the classical analytic signal 
to the complex case: it is also necessary to generalise the 
envelope and phase concepts, and this can only be done 
by interpreting the original and orthogonal complex signals 
as a pair. In this paper we have only one way to do this: 
by representing the pair of complex signals as a quaternion 
signal. This arises naturally from the third method above 
for creating an orthogonal signal, but also from the Cayley- 
Dickson construction of a quaternion as a complex number 
with complex real and imaginary parts (with different roots of 
— 1 used in each of the two levels of complex number). 

Extension of the analytic signal concept to 2D signals, that 
is images, with real, complex or quaternion-valued pixels is 
of interest, but outside the scope of this paper Some work has 
been done on this, notably by Biilow, Felsberg, Sommer and 
Hahn fSf-plj. The principal issue to be solved in the 2D case 
is to generalise the concept of a single-sided spectrum. Hahn 
considered a single quadrant or orthant spectrum, Sommer 
et al considered a spectrum with support limited to half the 
complex plane, not necessarily confined to two quadrants, but 
still with real sample or pixel values. 

Recently, Lilly and Olhede p2) have published a paper on 
bivariate analytic signal concepts without explicitly consid- 
ering the complex signal case which we cover here. Their 
approach is linked to a specific signal model, the modulated 
elliptical signal, which they illustrate with the example of a 
drifting oceanographic float. The approach taken in the present 
paper is more general and without reference to a specific signal 
model. 
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II. Construction of an orthogonal signal using 

CLASSICAL COMPLEX TRANSFORMS 

We first review and examine the notion of orthogonality 
in the case of complex signals before we consider how to 
construct a signal orthogonal to a given complex signal. 

Definition 1: Two complex signals x{t) and y{t) are or- 

thogonal if: / x{t)y*{t) At — 0, where * denotes the complex 

— oo 

conjugation. 

Orthogonality is invariant to multiplication of either signal 
by any real or complex constant. This is evident from the 
definition, since if the integral is zero for a given pair of x{t) 
and y{t) it must also be zero for C,y{t) or C^(^)' where C, 
is a complex constant. Therefore, given a pair of orthogonal 
signals, one of the signals may be rotated by an arbitrary angle 
in the complex plane (constant with respect to time), and still 
remain orthogonal to the other. Obviously the same applies 
for a real constant, so scaling one or both of the signals does 
not alter their orthogonality. 

Consider the integral in Definition [T] separated into real and 
imaginary parts; 



{Xr(t)yr(t)+X,{t)y^{t))^^t = Q (1) 

{x,(t)yr{t)-Xr{t)y^{t))dt^Q (2) 



where Xr{t) (resp. yr{t)) is the real part of x{t) (resp. y{t)) 
and Xi{t) (resp. yi{t)) is the imaginary part of x{t) (resp. 
y{t)). There are two significantly different ways in which these 
integrals can vanish. One way is for all four of the products 
within them to vanish when integrated: 



-oo 

oo 



Xr{t)yr{t) At^O J x,{t)y,lt) At =0 

— oo 
oo 

x,{t)yr{t) di = / Xr{t)y^{t) At =0 



This can only happen if both complex signals have an imag- 
inary part orthogonal to their real parfl In this case only, 
orthogonality is invariant to complex conjugation of either or 
both of the complex signals. 

In more general cases the integrals in Q and (|2]i vanish 
because the product terms within them integrate to equal 
values and appropriate signs (opposite signs in the first case, 
the same sign in the second). 

Of course, a hybrid of these possibilities is also possible, 
and indeed occurs in the next section: one integral vanishes 
because each of the product terms within it integrates to zero, 
and the other vanishes because the product terms integrate to 
the same value with appropriate sign, and cancel. 
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Fig. 1. Construction of the orthogonal signal using a pair of Hilbert 
transforms. 



A. Hilbert transform 

An orthogonal signal may be constructed using a pair of 
classical Hilbert transforms as shown in Figure [T] 

Lemma 1: Given two arbitrary real-valued functions f{t) 
and git) the following holds, where T-L [x{t)] is the Hilbert 
transform of x{t), provided the Hilbert transforms exist: 



f{t)n [g[t)] At+ H [f{t)] g{t) At = Q (3) 



Proof: We make use of the inner product form of Parse- 
val's Theorem y_3( §27.9 equation (27.34)] or | ,14, Theorem 
7.5]: 



x{t)y*{t) At 



X{v)Y*{v) Av 



(4) 



where X{v) is the Fourier transform of x{t) and similarly 
for y{t). Note that, if y{t) is real, as it is in what follows, 
Y{v) must still be conjugated. We also make use of the 
Fourier transform of a Hilbert transform | [T5 | (2.13-8, p68)]: 

F[l-i[x(t)]] = -is\gn{v)X{v), where X{u) = T[x{t)]. 
Using (|4]), we rewrite Q as: 

oo 

F*{u)isign{u)G{v) di^+y" jsign(zy)F(zy)G*(i/) dz^ = 

— OO 

which may be simpHfied to: 

oo 

i I sign{v) [F*{v)G{v) + F{v)G*{v)] Av = Q 



Since F*{v)G{v) and F{v)G*{v) are conjugates, their imagi- 
nary parts cancel and the term inside the brackets is real. Since 
F{v) and G{v) are the Fourier transforms of real signals, their 
real parts are even functions of v. Therefore when multiplied 
by the signum function, they integrate to zero. ■ 
Proposition 1: A complex signal o{t) constructed from 
an arbitrary complex signal z{t) using two classical Hilbert 
transforms operating independently on the real and imaginary 
parts, as represented in Figure [T] and below: 



o{t)=H[z,{t)]-iH [zm 
is orthogonal to z{t), that is: 



(5) 



'An example of such a signal can be created numerically by taking four 
singular vectors from a singular value decomposition of a real matrix, since 
the singular vectors are mutually orthogonal. 
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/ z{t)o* (t) dt ^ 

— oo 

Proof: First consider the real part of the integral: 



[Zr{t)Or{t) + Z,{t)o,{t)] dt = 



and write it as two integrals with the real and imaginary parts 
of o{t) replaced by their equivalents from (|5]|; 



Zrit)H [Zrit)] dt - Z,{t)H [Z.,{t)] dt = 



z{t)- 



Z{v) 



T- 



-o*{t) 



Both integrals evaluate to zero as a consequence of orthogo- 
nality between: Zr{t) and its Hilbert transform; and Zi{t) and 
its Hilbert transform. Now consider the imaginary part of the 
integral: 



[z,{t)Or{t) ~ Zr{t)0,(t)\ dt = 

and replace the real and imaginary parts of o[t) as before: 

oo 

{Zr{t)\ dt + J Zr{t)H [Z,{t)] dt = 
— oo 

This is true by Lemma [T| 

■ 

B. Complex Fourier transform pair 

An alternative method to construct an orthogonal signal 
is shown in Figure |2] using a complex Fourier transform 
pair operating directly on the complex signal, rather than 
Hilbert transforms operating on the real and imaginary parts 
separately. The Fourier spectrum of z{t) is multiplied by 
a signum function times i, equivalent to multiplication of 
positive frequencies by e™/"^ and negative frequencies by 
g-i7r/2 rpjjjg jj^g jjjg effect of rotating positive frequency 
coefficients by 7r/2 and negative frequency coefficients by 
— 7r/2. An inverse Fourier transform then yields the conjugate 
of the orthogonal signal. 

Proposition 2: A complex signal o(t) constructed from 
an arbitrary complex signal z{t) using a complex Fourier 
transform as represented in Figure [2] and below, is identical to 
the orthogonal signal in Proposition [T| equation (|5]l. 



o*{t)=T-^ [-isign{v)Z{v)] 



(6) 



where F denotes a Fourier transform and Z{v) = F [z{t)]. 

Proof: From equation o*{t) = U [zr{t)] + iU [zi{t)]. 
Taking the Fourier transforms of the Hilbert transforms (see 
Lemma [T]), we have: 

T[o*{t)\ — —isig-a{v)Zr{i') — isign(i^)Zi(i/) 
= —i sig'a{v)Z{v) 

from the linearity of the Fourier transform. 



'i sign(i') 



Fig. 2. Construction of tlie ortliogonal signal using a complex Fourier 
transform. 



C. Discussion 

We have shown in this section how to construct a complex 
signal o(t) orthogonal to an arbitrary complex signal z{t) 
using classical Hilbert or complex Fourier transforms. The 
notion of analytic signal, however, requires more than just 
an orthogonal signal — generalising from the classical case, 
we must treat the original and orthogonal complex signals 
as a pair in order to define the generalisations of envelope 
and phase. In the next section we introduce the quaternion 
Fourier transform, and show that this permits us to construct a 
hyperanalytic signal from which both the original signal and 
the orthogonal signal defined in the current section may be 
extracted. Using the quaternion hyperanalytic signal we then 
show how to construct the envelope and phase in the complex 



III. Construction of an orthogonal signal using a 

QUATERNION FOURIER TRANSFORM 

In this section, we wiU be concerned with the definition and 
properties of the quaternion Fourier transform (QFT) of com- 
plex valued signals. Before introducing the main definitions, 
we give some prerequisites on quaternion valued signals. 

A. Preliminary remarks 

Quaternions were discovered by Sir W.R. Hamilton in 1843 



1 16|. Quaternions are 4D hypercomplex numbers that form a 
noncommutative division ring denoted H. A quaternion q G H 
has a Cartesian form: a + bi + cj + dk, with a, 6, c, c? e M 
and i,j,k roots of —1 satisfying = = = ijk = 
~l. The scalar part of q is a: S{q) — a. The vector part 
of q is: V{q) = q — S{q). Quaternion multiplication is not 
commutative, so that in general qp ^ pq for p,q E H. The 
conjugate of q is q = a — bi — cj — dk. The norm of q is 
Ikll = = (a^+^^+c^+ci^) = gg. A quaternion with a = 
is called pure. If \q\ ~ 1, then q is called a unit quaternion. The 
inverse of q is q^^ = q/\\q\\- Pure unit quaternions are special 
quaternions, among which are i, j and k. Together with the 
identity of H, they form a quaternion basis: {l,i,j,k). In 
fact, given any two unit pure quaternions, fi and ^, which are 
orthogonal to each other (i.e. S{fj,^) = 0), then /i^) 
is a quaternion basis. 

Quaternions can also be viewed as complex numbers with 
complex components, i.e. one can write q = zi + jz2 in the 
basis (l,i,j,fc) with zi,Z2 G C*, i.e. z^ = ^(za) + i'^{za) 
for a = 1,2^ This is called the Cayley-Dickson form. 
Now, it is possible to define a generalized Cayley-Dickson 

-In the sequel, we note as a complex number with root of —1 being 
/i. Note that these are degenerate quaternions. 
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form of a given quaternion using an arbitrary quaternion 
basis {1, fi, ^, fi^). A quaternion q can then be written as 
q = Wi+ where wi,W2 G C^. 

Any quaternion q also has a unique polar Cayley-Dickson 
form flTl given by: q — Aexp{Bj), where A is the complex 
modulus of q and B is its complex phase. Expressions of A 
and B can be found in fn\ . This notation will be of use in 
Section |Vl 

B. ID Quaternion Fourier transform 

In this paper, we use a ID version of the (right) QFT first 
defined in discrete-time form in [18|. Thus, it is necessary to 
specify the axis (a pure unit quaternion) of the transform. So, 
we will denote [] a QFT with transformation axis /x. This 
will become clear from the definition below. We now present 
the definition and some properties of the transform used here. 

Definition 2: Given a complex valued signal z{t) — Zr{t) + 
$Zi{t), its quaternion Fourier transform with respect to axis /j, 
is: 



+ 00 



and the inverse transform is: 



+ 00 



(7) 



(8) 



where in both cases, the transform axis fi can be chosen to be 
any unit pure quaternion orthogonal to ^ so that (1, ^, /x, ^fi) 
is a quaternion basis. 

Property 1: Given a complex signal z{t) — Zr{t) + ^Zi{t) 
and its quaternion Fourier transform denoted Z^j,{v), then the 
following properties hold: 

• The even part of Zr{t) is in 3fi(Z^(z^)). 

• The odd part of Zr{t) is in 

• The even part of Zi{t) is in (;/)). 

• The odd part of Zi{t) is in 

Proof: Expand (|7]i into real and imaginary parts with 
respect to ^, and expand the quaternion exponential into cosine 
and sine components: 

+00 

Z^^{v) = j [zr{t) + [cos{2TTut) - /X sin(27ri^t)] At 



+ 00 



+ 00 



Zrit) cos{2TTvt)At — jj, I Zr{t) shi{2'n:vt)dt 



+00 



+00 

i{t) cos{2TTiyt)dt — ^fj, J Zi{t) svsi{2TTvt)&t 



from which the stated properties are evident. ■ 
These properties are central to the justification of the use 
of the QFT to analyze a complex valued signal carrying com- 
plementary but different information in its real and imaginary 
parts. Using the QFT, it is possible to have the odd and even 
parts of the real and imaginary parts of the signal in four 



different components in the transform domain. This idea was 
also the initial motivation of Biilow, Sommer and Felsberg 
when they developed the monogenic signal for images [SI- 
IT^. Note that the use of hypercomplex Fourier transforms 
was originally introduced in 2D Nuclear Magnetic Resonance 
image analysis |19J , | ,20J . 

We now turn to the link between a complex signal and the 
quaternion signal that can be uniquely associated to it. 

Property 2: A one-sided QFT X{v) with X(v) ^ 
(i;eSlv < 0) corresponds to a quaternion valued signal in the 
time domain. 

Proof: A complex signal with values in has a QFT^ 
with the symmetry properties listed in property [T] Now, 
cancelling its negative frequencies is equivalent to adding a 
signal with a QFT with the following properties: part is 
odd, 5^ is odd, 3^ is even and is even. By linearity of 
the inverse QFT, the former two parts of the QFT produce a 
complex signal in the time domain with 3? and 3^ parts, while 
the two latter produce a quaternion signal in the time domain 
with 3^ and parts. Together, and again by linearity, the 
resulting signal in the time domain is a full 4D quaternion 
signal. ■ 
Property 3: Given a complex signal x{t), one can associate 
to it a unique canonical pair corresponding to a modulus and 
phase. These modulus and phase are uniquely defined through 
the hyperanalytic signal, which is quaternion valued. 

Proof Cancelling the negative frequencies of the QFT 
leads to a quaternion signal in the time domain. Then, any 
quaternion signal has a modulus and phase defined using its 
CD polar form. ■ 

C. Convolution 

We consider the special case of convolution of a complex 
signal by a real signal. Consider g and / such that: g : M+ ^• 
and / : M+ ^ M. Now, consider the QFT^ of their 
convolution: 

[g * fit)] ^ J J sir) fit - r)dre-2^-"^Mi 



— oo — oo 

+ 00 +00 



giT)e-^'^'"''-''+^^fit')dTdt' 



— OO — OO 
+ 00 



+ 00 



(^)g-2M^..rj^ / fit')e-^'''"'*'dt' 



= 9 



= T^[git)]T^[fit)] 

(9) 

Thus, the definition used for the QFT here verifies the convo- 
lution theorem in the considered case. This specific case will 
be of use in our definition of the hyperanalytic signal. 

D. The quaternion Fourier transform of the Hilbert transform 

It is straightforward to verify that the quaternion Fourier 
transform of a real signal xit) — ^ is — sign(i^), where /x 
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is the axis of the transform. Substituting x{t) into (|7]), we get: 

+00 



' 1 ■ 






_7rt 







t 



At 



and this is clearly isomorphic to the classical complex case. 
The solution in the classical case is — isign(z^), and hence in 
the quaternion case must be as stated above. 

It is also straightforward to see that, given an arbitrary real 
signal y{t), subject only to the constraint that its classical 
Hilbert transform T-L [y{t)] exists, then one can easily show that 
the classical Hilbert transform of the signal may be obtained 
using a quaternion Fourier transform as follows: 



H[y{t)]=T^^[-~^ls\gn{v)Y^{v)] 



(10) 



where Y^{v) — Ffj, [y{t)\. This result follows from the 
isomorphism between the quaternion and complex Fourier 
transforms when operating on a real signal, and it may be 
seen to be the result of a convolution between the signal y{t) 
and the quaternion Fourier transform of x{t) — Note that 
/i, and Y^{v) commute as a consequence of y{t) being real. 

IV. The hyperanalytic signal as a quaternion 

SIGNAL WITH A ONE-SIDED SPECTRUM 

We define the hyperanalytic signal h{t) by a similar ap- 
proach to that originally developed by Ville |TJ. The following 
definitions give the details of the construction of this signal. 
Note that the signal z{t) is considered to be non-analytic in 
the classical (complex) sense, that is its real and imaginary 
parts are not orthogonal. However, the following definitions 
are valid if z(t) is analytic, as it can be considered as a 
degenate case of the more general non-analytic case. 

Definition 3: Consider a complex signal z{t) — Zr{t) + 
£,Zi{t) and its quaternion Fourier transform Zn{i^) as defined 
in Definition |2] Then, the hypercomplex analogue of the 
Hilbert transform of z{t), is as follows: 



(11) 



where the Hilbert transform is thought of as: H/j, [z{t)] = 
p.v. (z * where the principal value (p. v.) is understood in 
its classical way. This result follows from equation ( [TO] i and 
the linearity of the quaternion Fourier transform. Notice that 
the result will take the same form as z(t), namely a quaternion 
signal Or{t)+^Oi{t) isomorphic to the complex signal o*(t) = 
0r{t) + iOi{t). To extract the imaginary part, the vector part of 
the quaternion signal must be multiplied by — ^. An alternative 
is to take the scalar or inner product of the vector part with ^. 
Note that /j, and {v) anticommute because fi is orthogonal 
to ^, the axis of Z^{v). Therefore the ordering is not arbitrary, 
but changing it simply conjugates the result. 

Definition 4: Given a complex valued signal z{t) that 
can be expressed in the form of a quaternion as z{t) = 
Zr{t) + ^Zi{t), and given a pure unit quaternion /j, such 
that (1, ^, fi, ^/i.) is a quaternion basis, then the hyperanalytic 
signal of z{t), denoted h{t) is given by: 



where "H^ [z{t)] is the hypercomplex analogue of the Hilbert 
transform of z{t) defined in the preceding definition. The 
quaternion Fourier transform of this hyperanalytic signal is 
thus: 

H^{i^) = Z^{u) - sign{v)Z^{v) 
= [1 + sign(i/)] 

which is a direct extension of the 'classical' analytic signal. 

This result is unique to the quaternion Fourier transform 
representation of the hyperanalytic signal — the hyperan- 
alytic signal has a one-sided quaternion Fourier spectrum. 
This means that the hyperanalytic signal may be constructed 
from a complex signal z{t) in exactly the same way that 
an analytic signal may be constructed from a real signal 
x{t), by suppression of negative frequencies in the Fourier 
spectrum. The only difference is that in the hyperanalytic case, 
a quaternion rather than a complex Fourier transform must be 
used, and of course the complex signal must be put in the form 
z{t) = Zr{t) + $Zi{t) which, although a quaternion signal, is 
isomorphic to the original complex signal. 

A second important property of the hyperanalytic signal is 
that it maintains a separation between the different even and 
odd parts of the original signal. 

Property 4: The original signal z{t) is the simplex part ||2T| 
Theorem 1], |j22j § 13.1.3] of its corresponding hyperanalytic 
signal. The perplex part is the orthogonal or 'quadrature' 
component, o{t). They are obtained by: 



z{t) = \ (hit) - ^hm) 
o{t)^^-{h{t)+^hm 



(13) 
(14) 



Proof: This follows from equation ( [T2] l. Writing this in 
full by substituting the orthogonal signal for T-L^ [z{t)]: 

h{t) = z{t) + fl0{t) ^ Zr{t) + iz,{t) + flOr{t) + ^1^0^it) 

and substituting this into equation ( [T3| ), we get: 

1 / Zr{t)+ ^Z, (t) + flOr (t) + fl^O, (t) \ 

2 V [Zr{t) + ^Z,{t) + HOr{t) + ^l^.O^{t)] ^ ) 



Z{t) 




h{t) = z{t) + fiH^j, [z{t)] 



(12) 



and since ^ and /x are orthogonal unit pure quaternions, ^/x = 

Zr{t) + iZi{t) + fJ,Or{t) + fl^Oi{t)\ 
Zr{t) + iZi{t) - flOr{t) ~ fl^O,{t) J 

from which the first part of the result follows. Equation 
( [l4l ) differs only in the sign of the second term, and it is 
straightforward to see that if h{t) is substituted, z{t) cancels 
out, leaving o{t). ■ 

V. Phase and envelope from quaternion 

REPRESENTATION 

As was shown in § II-A and Figure [T] it is possible to 
construct an orthogonal signal using classical Hilbert trans- 
forms operating independently on the real and imaginary parts 
of a complex signal z{t). However, to extend the analytic 
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signal concept to the hyperanalytic case, we need a method 
to construct the complex envelope (also known as the in- 
stantaneous amplitude) and the complex phase (also known 
as the instantaneous phase) from a pair of complex signals 
(the original complex signal z{t), and the complex orthogonal 
signal o{t)). 

A complex envelope may be obtained, apart from the signs 
of the real and imaginary parts, and therefore the quadrant, 
from the moduli of the analytic signals of the real and 
imaginary parts of the original signal. The resulting envelope 
will always be in the first quadrant and therefore in general it 
will not be a correct instantaneous amplitude. 

In this paper the method we present is based on quaternion 
algebra which is a natural choice given that a quaternion 
may be expressed as a pair of complex numbers in the so- 
called Cayley-Dickson form, and, as we have shown earlier, 
the hyperanalytic signal may be computed using a discrete 
quaternion Fourier transform. 

A. Complex envelope 

In a previous paper fSl we showed how the envelope of the 
hyperanalytic signal could be obtained from a biquaternion 
Fourier transform, using simply the complex modulus of 
the biquaternion hyperanalytic signal. This approach has a 
limitation, namely that the envelope is always in the right 
half of the complex plane. (This is because the modulus of 
the biquaternion hyperanalytic signal is obtained as the square 
root of the sum of the squares of the components, and the 
complex square root is conventionally in the right half-plane.) 
A second, more serious, limitation of the biquaternion Fourier 
transform is that under some circumstances the complex 
envelope can vanish because the biquaternion hyperanalytic 
signal is nilpotent [23]: that is it has values with zero modulus, 
even though the biquaternion values are non-zero. 

The mathematical key to constructing the complex envelope 
from a quaternion hyperanalytic signal is given in a paper 
written in 2008 |17J . In that paper we set out a new polar 
representation for quaternions that has a complex modulus and 
a complex argument: 

g = ^exp(Bj) (15) 

where A and B are complex (more accurately they are 
degenerate quaternions of the form it; + + Oj + Ofc). 

The representation as published is based on the standard 
quaternion basis (1, i, j, fe), although clearly it could use any 
other basis. For simplicity of numerical computation, we use 
the standard basis. 

This representation of a quaternion is not trivial to compute 
numerically and one of the main problems we had in arriving 
at the results presented in this paper was to refine the com- 
putation of the representation in ([TSj to minimise numerical 
problems. We do not discuss these issues here, as they are 
solely concerned with correct numerical implementation of 
([TSj. The code is open source, and the numerical issues are 
documented in the code \2A^ Function: cdpolar . m]. It is 
also necessary to use phase unwrapping on B to obtain a 
continuous envelope and phase without discontinuities when 
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Fig. 3. Classical envelope of an amplitude modulated signal obtained as 
the modulus of the analytic signal. Blue: original modulated signal; green: 
orthogonal signal; red: envelope. 

computing A{t) and B{t) using ([TSj from a hyperanalytic 
signal h{t), due to an inherent ambiguity of sign in A 
and B |17, §2.1]. (In a similar way the complex envelope 
obtained from a biquaternion hyperanalytic signal may be 
phase unwrapped to obtain an envelope in all four quadrants.) 

Definition 5: Given a complex signal z{t), and its hyper- 
analytic signal h{t) constructed according to Definition |4] 
expressed in the standard quaternion basis {l,i,j,k), the 
instantaneous amplitude, or complex envelope, e(t), and the 
instantaneous phase, (j){t), are given by the polar form of a 
quaternion defined in ([T5|: 

h{t) = z{t) + o{t)j - e{t) exp{cj){t)j) 

where z{t), o{t), e{t), and 4){t) are complex. 

B. Complex modulation of a real carrier 

A classic result in the case of the analytic signal concerns 
amplitude modulation, where the envelope of the analytic 
signal (that is its modulus) allows the modulating signal to 
be recovered, subject to certain constraints such as the depth 
of modulation |T, §7.16]. Figure |3] shows an example. In this 
section we extend this result to the complex case. 

Property 5: Given a signal z{t) consisting of a real si- 
nusoidal 'carrier' f{t) = A cos(27rti;/; + 0), and a complex 
'modulating' signal g{t): 

z{t) - g{t) fit) 

where the frequency content of g{t) is below the frequency 
of the carrier, lo, then the instantaneous phase 0(t) as defined 
in Definition |5] is real and equal to the instantaneous phase 
of f{t), that is 27ra;t + 0, and the envelope e(t) is a scaled 
version of the complex modulating signal, Ag{t). 

Proof Separate z{t) into real and imaginary parts: 

Zr{t) = gr{t)Acos{2TrLut + 9) 
z^{t) = gi{t)Acos{2nujt + 9) 
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By Bedrosian's Theorem |]25| Q §7.12 and Table 7.7.1], the 
Hilbert transforms of the real and imaginary parts of z{t) are: 

n [zr{t)] = gr{t)Asm{2TTiot + 0) 
n[z,{t)] =5i(t)Asin(27rw< + 61) 

We can therefore write the hyperanalytic signal as follows: 

h{t) = Zr{t) + iZi{t)+jn [Zr{t)] + kn 

= g{t)A[cos{2Trujt + 6) + j sin(27rwt + 9)] 
= g{t)Aeiip{j{2TTUJt + 0)) 

and comparison with Definition |5] shows that e(t) = Ag{t) 
and (p{t) = 2Trojt + 6. ■ 
Note that we have considered above only the case of a real 
sinusoidal carrier modulated by a complex modulating signal, 
and that scaling of the carrier cannot be separated from scaling 
of the modulating signal. A more general case would be a 
complex carrier, for example, a complex exponential, or a 
single sinusoid rotated in the complex plane, modulated by 
a real or complex modulating signal, but we now show that 
this is ambiguous - the complex nature of the carrier cannot be 
separated from the complex nature of the modulating signal. 
Consider a real sinusoidal carrier f{t) modulated by a complex 
modulating signal g(t), yielding a modulated signal m{t): 

m{t)^f{t)g{t) 

Now rotate the carrier by an angle ^ in the complex plane so 
that it oscillates along a line in the complex plane at angle tp 
to the real axis, instead of along the real axis itself: 

m{t) = [/(t) exp(zV')] <?W 



Clearly, this modulated signal is indistinguishable from the 
following, in which the complex modulating signal g{t) is 
instead rotated through the angle ip: 

m{t) = f{t) [g{t) exp(iV)] 

We conclude that the hyperanalytic signal cannot distinguish 
between a real carrier modulated by a complex signal and a 
complex carrier modulated by a (different) complex signal: 
it suffices to consider the former case, and assume that the 
modulating signal only is complex. 

VI. Properness 

In this section, we study the properness properties of the 
hyperanalytic signal. We make use of the term proper here 
for deterministic signals, by extension of the concept used 
for quaternion random functions in | |26[ . The term proper 
is also widely used for complex random signals when the 
real and imaginary parts are orthogonal and have the same 
amplitude. Our definition of properness is based on p7) . 
We use the properties of the correlation matrix in the real 
(Cartesian) representation of quaternions. As explained in p7) , 
a quaternion can be represented as a four dimensional vector 
q — [a,b,c^dY' ■ Properness can be observed through the 
pattern of the covariance matrix Tq that contains the cross- 
correlations between the four components of the quaternion. 



Here, we are interested in the covariance (correlation at zero- 
lag) of the hyperanalytic signal h{t). Assuming that the 
original signal z{t) takes values in C^, then h{t) (obtained 
using a quaternion Fourier transform with axis /x) can be 
written as: 

h{t) = Z(t) + o{t)n = Zr[t) + $,Z,{t) + {Or{t) + |o,(t)) /X 

The covariance matrix of h{t), denoted Th, is thus given as: 





Cz^Zi 


Cz^Or- 


Cz^o 


CziZr 


^ZiZi 


CziOr 


CziO 


Co^Zr 


Co^ Zi 




Co^o 


C-OiZr 


^OiZi 


CoiOr 


CoiO 



Now, given a complex signal z{t) taking values in and 
its hyperanalytic (quaternion valued) signal h{t), then the 
following stands. 

Property 6: If z{t) with values in is an improper 
complex signal, then its hyperanalytic signal h{t), computed 
using a quaternion Fourier transform with axis /x is C- 
proper, providing that (1,^, is a quaternion basis. As 
a consequence, the covariance of h{t) is: 

a 7 (3' 

7 w /3 

0/3 a -7 

/3 -7 w 



with: 



a 



z^{t)dt 



ol{t)dt 



zt{t)dt= / ofit)dt 



Zi{t)Or{t)dt = / Zr{t)oi{t)dt 



'~f ~ J Zr{t)zi{t)dt = — J 0r{t)0i{t)dt 

— oo — oo 

Proof: We need to prove the zeroing of the terms of F/j 
as well as equalities between some of its terms. Obviously 
Tfi is a real symmetric matrix, so we only need proofs for 
the diagonal and upper-diagonal elements. We start with the 
diagonal elements. 

First, recall that, due to Parseval's theorem: 



Zritfdt 



\Zr{v)\^dv 



Now by definition of the orthogonal signal o{t), the Fourier 
transform of its real part is: 

Or{y) = T K(i)] ^T{n {z.,{t)\\ = -i^Y^^{v)Z,{v) 

It follows directly that |Or(i/)p = and using Parse- 

val's theorem (which is also valid for Orit)), the following is 
true: 



zl{t)dt= / ol(t)dt^ \ U{Zr(t)fdt 
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A similar reasoning is valid to show that: 

/oo noo 
z,{tfdt= 1 ol{t)At 
-oo J —oo 

with the small difference that Zi{t) — [zi{t)], but the sign 
has no effect on the reasoning due to the modulus. 

We now turn to off-diagonal elements. As z{t) is an 
improper signal, its real and imaginary parts are not orthogonal 
to each other, so the following is true: 



We now show that in the case where z{t) is improper, then 
/3 9^ 0. We start by recalling that: 



Z,{t)Or{t)dt = / Zi{v)0*{u)dv ^ / Or{v)Z*{v)dv 



— OO 

oo 



Zr{t)o,{t)dt ^ / Z^(j/)0.*(t/)di/ = / Zriiy)0*{iy)di^ 



7 = y Zr{t)z,{t)dt^O 

— oo 

Now, Parseval's theorem allows us to write: 

oo 

Zr{t){z*{t)ydt = J Zr{v)Z,{-v)du 
— oo 

where Zr{v) and Zi{v) are the classical Fourier transforms of 
Zr{t) and Zi{t). We also have that: 

o^{t) = -niz.it)] 

which implies that: 

Or{v) = -i%\'gi\{y)Zr{y) 
0,(y)^ isx%x,(v)Z,{v) 

Therefore, we can write the following equaUty: 



Now, using the expressions for Or(y) and Oi(y) given above, 
we have that: 



z,(t)or{t)dt = i \ iiYgx^{y)Z,{v)Z*\v)dv 



oo oo 
J Zr{t)o^{t)dt = j S,lg-a{v)Zr{y)Z*{v)dv 

— OO — oo 

But, remember that Zr{t), Or{t), Oi{t) and Zi{t) are real valued, 

oo oo 

and thus / z^{t)Or[t)dt G M and / Zr{t)oi{t)dt G M. As a 

— oo — oo 

consequence, the two integrals equal to p may be written: 



Zi{t)Or{t)dt = — 



-oo 
oo 



Zr{t)o,{t)dt ^ 3 



sigii{v)Z,{y)Z*{v)dy 



sign{v)Zr{v)Z*^{v)dv 



Or{v)Oi{-v)dv 

j -isign{v)Zr{v)is\gn{-v)Z,{-u)dv 



As these integrals are all equal, we simply need to show that 
one of them is not equal to zero for our purpose. Now, let us 
consider the following quantity: 



W{y) = S 



Zr{v)Zi{ — v)dv = — J Zr{t)zi{t)dt 

— oo 

This shows, due to Parseval's theorem, that: 

oo 

Or{t)oi{t)dt — — Zr{t)zi{t)dt = —J 



Zr{y)Z*[v)dv 



One can see by direct calculation (using Hermitian symmetry 
of the classical Fourier transform) that: 



W{-v) = 3? 



Then, we have already shown in the proof of Proposition [T] 
that: 

oo oo 

j Zr{t)or{t)dt = and j z,{t)o,{t)dt = 



due to orthogonality between each signal and its Hilbert 
transform. 

Now, it was also shown in Proposition [T] that, due to 
orthogonality between z{t) and o{t) the following stands: 



Z,{t)Orit)dt = / Zrit)o,{t)dt ^ l3 



Zr{-v)Z*{-v)dv 

Z*{y)Z,{v)dy 

Zr{v)Z*{v)dv 



which shows that W{v) is odd. Now, as a consequence, the 
following is true: 



sign{v)Zr{v)Z*{v)diy 
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And this makes it impossible to cancel out, except in the very 
special case where Zr{i^) and/or Zi{i') are equal to zero for 
any v. Thus the following is true: 



sigTi{iy)Zr{iy)Z*{u)du 



Zr{t)o^{t)dt = /3 7^ 



This completes the proof. ■ 

In the most 'classical' case {z{t) taking values in C*), this 
means that if one computes the corresponding hyperanalytic 
signal using a quaternion Fourier transform with axis j, then 
the hyperanalytic signal h{t) will be -proper 

A way to figure out the meaning of the C-properness of the 
hyperanalytic signal is to see it as a pair of complex signals, 
both improper and orthogonal to each other (in the sense of 
the complex scalar product) but not jointly proper 

We now consider the case when the original signal is already 
a complex analytic signal. This is a degenerate case, and the 
consequences for properness properties are given as follows. 

Property 7: If z{t) with values in is a proper complex 
signal, then its hyperanalytic signal h{t), computed using a 
quaternion Fourier transform with axis fi is C^-superprope:]^ 
provided that (1.^, is a quaternion basis. As a conse- 
quence, the covariance of h{t) reads: 



a 








/? 





a 










fi 


a 





/3 








a 



zf.{t)dt= / ol{t)dt^ / zf{t)dt= / o^{t)dt 



with: 



and (3 being as defined in Property |6] 

Proof: Referring to Property [6] we need to prove that: 
a — u! and 7 = 0. 

First, as z{t) is a complex analytic signal, it is proper and 
so: 



Zr{t)z,{t)dt = 7 = 



Then, it is also a well-known property of analytic signals that 
their real and imaginary parts have the same magnitude. Here, 
this means that: 



z^(t)dt 



zfit)dt 



and consequently a ^ uj. This completes the proof. 

Note that we have only highlighted here the properness 
properties of the hyperanalytic signal and left for future work 
the possiblity of taking advantage of these properties in its 
processing. ■ 

^^We use the term C^-superproper for a C^-proper quaternion valued signal 
obtained from a complex proper signal, i.e. for which 7 = in the covariance 
matrix. 




Fig. 4. Sinusoidal carrier modulated by a complex helix and two Gaussian 
pulses. Blue: original signal; green: orthogonal signal; red: complex envelope. 



VII. Examples 

We present here an example signal, with its complex en- 
velope and phase computed using a hyperanalytic signal. The 
hyperanalytic signal is computed from a one-sided quaternion 
Fourier transform, with the negative frequencies suppressed 
using the discrete-time scheme first published by Marple p8) . 

Figure |4] shows a complex analogue of the signal shown 
in Figure |3] The sinusoidal carrier is amplitude modulated 
by two Gaussian pulses, and modulated in angle (in the 
complex plane) by multiplication with a complex exponential. 
The figure shows an orthogonal signal and complex envelope, 
computed from a hyperanalytic signal constructed using a 
quaternion Fourier transform, with the envelope extracted as 
described in Definition [5] Figure [5] shows the real and imagi- 
nary components of the signals shown in Figure |4] as well as 
the original sinusoidal carrier and the extracted complex phase 
(j){t) (in this case the imaginary part is negligible, as expected 
from Property |5] and is not shown). The envelope can be seen 
to pass through all four quadrants. Phase unwrapping was used 
to eliminate discontinuities in the complex envelope, operating 
on the angle of the complex envelope and constructing the final 
result from the modulus of the complex envelope combined 
with the phase-unwrapped angle. The QTFM toolbox | |24| 
was used for the computations, and in particular, the function 
cdpolar.m was used to extract the complex envelope from 
the hyperanalytic signal. 

VIII. Conclusions 

We have shown in this paper how the classical analytic 
signal concept may be extended to the case of the hyperan- 
alytic signal of an original complex signal, constructing an 
orthogonal complex signal by one of three equivalent methods. 
The quaternion method yields an interpretation of the hyper- 
analytic signal as a quaternion signal which leads naturally to 
the definition of the complex envelope. We have shown that 
the quaternion polar representation developed by the authors 
in 2008 permits us to recover the phase of a sinusoidal carrier 
modulated by a complex signal with frequency content lower 
than the carrier frequency, extending the classical result of 
amplitude modulation. 
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Fig. 5. Top: real parts of signals in Figure |4] Middle: imaginary parts 
of signals in Figure |4] Lower: sinusoidal carrier (blue) and recovered canier 
phase (f>{t) (green), using Definition |5] 

We have presented the properness of the hyperanalytic 
signal, showing that for a given improper complex signal, its 
hyperanalytic signal is C-proper. 

As we stated earlier in the paper, our aim in this work 
has been to make some progress towards an analytic signal 
concept applicable to vector or quaternion signals. We believe 
the present paper indicates some of the problems to be 
overcome in doing this, including the need to work in a higher 
dimensional algebra than the quaternions, just as the results 
in this paper are most simply expressed using quaternion 
concepts rather than the concept of complex pairs. 

There is also scope to combine the ideas presented here 
with the ideas of monogenic signals (that is analytic signals 
in two dimensions, or the analytic signal of images) pO) . 
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